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Abstract 

Explicit solution of an infinite horizon optimal stopping problem for 
a Levy processes with a polynomial reward function is given, in terms of 
the overall supremum of the process, when the solution of the problem is 
one-sided. The results are obtained via the generalization of known results 
about the averaging function associated with the problem. This averaging 
function can be directly computed in case of polynomial rewards. To 
illustrate this result, examples for general quadratic and cubic polynomials 
are discussed in case the process is Brownian motion, and the optimal 
stopping problem for a quartic polynomial and a Kou’s process is solved. 
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1 Introduction 


Since the seminal work of Darling, Liggett and Taylor (1972), giving the solution 
to the optimal stopping problem for random walks, and reward functions of the 
form g(x ) = x + and g{x) = ( e x — 1) + , in terms of the distribution of the 
maximum of the random walk, it became clear the possibility of linking these 
two relevant problems in probability theory: the optimal stopping problem and 
the computation of the distribution of the overall maximum of a random walk. 
The natural question that this work posed was the possibility of extending 
these results to more general classes of processes, and to more general reward 
functions. 

The first results for Levy processes were obtained by Mordecki (21)l) 2a|b i. 
where the similar corresponding problems for arbitrary Levy processes are solved, 


*Universidad de la Republica, Facultad de Ciencias, Centro de Matematica. Igua 4225, 
11400 Montevideo, Uruguay, e-mail: mordecki@cmat.edu.uy 

'National Taras Shevchenko University of Kyiv, Department of Probability Theory, 
Statistics and Actuarial Mathematics. Volodymirska 60, 01601 Kyiv, Ukraine, e-mail: 
myus@univ. kie v. ua. 


1 









based on a discretization approximation argument, for the same reward func¬ 
tions, with the novelty of the consideration of the decreasing put reward g(x) = 
(K — e x ) + , that has a solution in terms of the overall infimum of the process. 
The first results for general payoffs were obtained by [Boyarchcnko and Leven-| 


dorskii (2002). Namely, using the technique of the Pseudo-Differential opera¬ 


tors, these authors obtained solutions to optimal stopping problems considering 
a large class of reward functions, making clear that the obtained previously re¬ 
sults were not based on particular properties of the payoff function, but only 
on the properties of the Levy processes. Their approach is analytic, based on 
the decomposition of an operator, that is in certain sense equivalent to the 
Wiener-Hopf factorization, and imposes certain restrictions on the class of Levy 
processes to which the results can be applied. For a general exposition of these 


results see also Boyarchenko and Levendorskii (2002a). Afterwards, Novikov 


and Shiryaev (2004) solved the optimal stopping problem for arbitrary random 
walks and reward functions of the form g(x) = (s + )", in terms of the Appell 
polynomials, and Novikov and Shiryaev] (2007) gave the solution to the problem 
with a power function reward with real and positive exponent, for both random 


walks and Levy process. Salminen (2007) applies the representation method for 


this problem (initiated in Salminen (1985)) finding the representing measure of 
the value function. More recently, Mishura and Tomashyk (20111 considered 
the optimal stopping problem for a general polynomial reward and a random 


walk. Alili and Kyprianou 

(20051 and 

Kyprianou and Surya 

(2005) obtained a 

new proof of the main results in 

Mordecki 

(2002a 

) and a generalization of the 

results in 

Novikov and Shiryaev 

(2004 

) for Levy process respectively, in both 


cases based on the strong Markov property of Levy processes. These contribu¬ 
tions were summarized in the monograph by Kyprianou (2006). On the way 
to the consideration of more general processes, Mordecki and Salminen ( 2007| ) 
obtained a representation of the value function for Hunt processes, that in the 
case of Levy processes give a representation in terms of the maximum of the 
process, and Christensen et al. (2013) exploited the excessive property of the 
maximum of a Markov process to obtain a verification theorem. It became then 
clear that the results were based on the probabilistic properties which random 
walks and Levy processes share, i.e. the independence and homogeneity of in¬ 
crements, and not on the particular form of the reward functions. Nevertheless, 
some particular reward functions admitted solutions in closed form. 

The approach that we use in this paper is the averaging problem , that was 
introduced in Surya (2007) (see also Surya (2007a)). The objective of the present 
paper is then twofold. We first present a theorem that summarizes and slightly 
improves the results of Surya (2007) and Christensen et al. (2013) in the case of 
Levy proce sses. The imp rovement consists in th e observation that t he ave raging 
function in 


Surya (2007) (or the function / in 


Christensen et al. 


(2013)) need 

not to be defined in the whole line, consequently the condition of this function to 
be negative on a certain set is not necessary. This allows to apply the result to 
larger classes of payoffs functions, what can be verified for certain polynomial 
rewards (see Remark [I]). The second objective of the paper is to apply the 
previous results to the class of general polynomial rewards. The main result 
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there is a simple algorithm to compute the averaging polynomial P of a given 
polynomial p. 

The content of the paper is as follows. In Section 2 we formulate the problem 
and prove the main results. In Section 3 we specializes to polynomial rewards. 
In Section 4 we present some examples: we discuss in detail the optimal stopping 
problem for Brownian motion and general quadratic and cubic polynomials, and 
also solve explicitely the optimal stopping problem for a quartic polynomial for 
a Kou process. 


2 Formulation of the problem and main results 

Let X = {X t } t > 0 be a Levy process defined on a stochastic basis B = (fi, P, F = 
(Pt)t>o, Pa) departing from Xq = x. For z G iR, the Levy-Khintchine formula 
states E 0 e zX * = ewith 

2 p 

tp(z) =az+ °—z 2 + / ( e zy - 1 - zh{y)) II {dy), (1) 

1 J R 

where a £ R, a > 0 and II (dy) that satisfies f R (l A y 2 )JJ(dy ) < oo conform 
the characteristic triplet (a, a, II) of the process. Here h(y) = is a 

truncation function. Given the stochastic basis B the set of stopping times is 
the set of random variables 


M = {t : H — > [0, oo] such that {r < t} G P t for all t > 0}. 


Observe that we allow the possibility r = oo, as for several optimal stopping 
problems, the optimal stopping time is within this class. A key role in the 
solution of the problem is played by the overall maximum of the process, defined, 
for r > 0 by 

M = sup{dG : 0 < t < e(r)}, 

where e(r) is an exponential random variable of parameter r > 0, and we assume 
e(0) = oo. We further assume thorough the paper that M is a proper random 
variable. This entails either that r > 0 or that X = {X t }t >o drifts to —oo 
when r = 0, and that P X (M > x) > 0, excluding the case of the negative of a 
subordinator, that gives M = x a.s. 

Given a non-negative payoff function g{x ), a process {AT ( }t>o departing from 
Xq = x adapted to a filtration F, and a discount factor r > 0, the optimal 
stopping problem consists in finding the value function V(x) and the optimal 
stopping rule r* such that 

V(x) = sup E x e~ rT g(X T ) = E x e~ rT *g(X T ,). (2) 

t£M 


Following Shiryaev (2008) we assume that the payoff received in the set {w: r (w) = 
oo} is 

limsupe -rt (/(A' t ). 
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In the present paper we are interested in problems with one-sided solution, i.e. 
such that the optimal stopping rule is of the form 


t* = inf{f > 0: X t > x*}, 


(3) 
g{x) = 


for some critical threshold x*. For this reason we assume that lim sup^^ 
lim^oo g(x) = 0. 

The averaging problem for optimal stopping, introduced by Surya| (2007) 
consists in finding an auxiliary function Q such that 


E X Q(M) = g(x) for all x G 


(4) 


where g is the payoff function of the problem and M the overall maximum. 
This approach, combined with the strong Markov property and invariance of 
increments of Levy process gives a fluctuation identity that allows to write the 
value function of the problem <[2j) in terms of M (see (10) in Lemma 1 below). 
Here we present a generalization of the results in Surya (2007). 


Theorem 1. Consider a Levy process X = {X t } t >o, a discount rate r > 0, and 
a reward function g: R. —» [0,oo) such that linr^-oo g(x) = 0. Assume that 


there exists a point x* and a non-decreasing function G* 
that 

E X G*(M) = g(x), for all x > x*. 


[x*, 00 ) 


such 


Define the function 


G{x) = 


G*(x), if x > x* 
0, if x < x* 


and the function V: R —> R. by 


If the condition 


V(x) = E X G(M), for all x G 


V(x) > g{x), for all x < x* 


(5) 


( 6 ) 


(7) 


is satisfied, then the optimal stopping problem ([2]) has value function V (x) in 
© , and © is an optimal stopping time for the problem. 


Remark 1. Compared to Theorem 5.3.1. in Surya (2007a), Theorem^above 


does not require the solution of the averaging problem for g and M to be found 
on the whole real line, but only on a certain set of the form [x*,oo). The relevant 
new condition to be verified on the set (— oo,x*) is ([7]). If the averaging function 
Q (satisfying Q) can be defined in the whole real line and it satisfies Q{x) < 0 
on the set (—oo,x*), then, condition Q follows (see Corollary [7]). Our function 
G is simply defined to be zero on this set. In Example \ f.l\ when a = — 1 we 
observe that Q is verified while the averaging function corresponding to Q 
takes positive values (for instance, Piifl) = see Figure 4.1). Furthermore, 


condition ([7|) is slightly more general than condition (b)(ii) in Theorem 2-4 in 
Christensen et al. (2013). 
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As usual in optimal stopping proofs we have to verify two statements: 


V(x) = E x e~ rT *g(X T ,) = E x e~ rT * g(X T ,). (8) 

V(x)>-E x e~ rT g(X T ) = -E x e- rT g(X T ), Vr e M. (9) 


These two statements are proved based on the following two lemmas which 


proofs follow essentially the respective proofs of Surya (2007) and Christensen 


et al. (2013) with the minor necessary modifications. 


Lemma 1. Consider a Levy process X, a discount rate r > 0, a reward function 
g, a threshold x*, an averaging function G*, and the extended function G, all 
this elements as in Theorem [7J Then, for any a > x* and x £ M, 


E,G(M)l {M > a} = E*e- rT “ 5 (X r Jl {Ta<oo} . (10) 

In particular, when a = x*, for t* in ®> and r > 0, we have 
E x G(M) 1 {m > x * } =E x e~ rT g(X T *) l{ r » <00 }. 
meanwhile, when r = 0, we have 

E x G(M) 1 {m > x , } = E x g(X T * )1{ T * <00 }. 

Proof. Consider, for a > x* , a hitting time of the form 

r a = inf{i > 0: X t > a}. 

As Levy processes satisfy the homogeneity property of increments in time and 
space, conditionally to the cr-algebra J> a , and on the set {r a < oo}, the process 
X s = X Ta+s — X Ta is independent of T Ta and has the same distribution as X (see 
Theorem 7, Chapter 4 in Skorokhod (1991)). We then consider two independent 
Levy processes X and X defined on a product probability space PxP. Consider 
first the case r = 0. We have 


E 


,G(M)l{M>a} — ExG( sup Xt) l{ Ta<00 } 

'■ 0<£<oo ' 

■G(x Ta + sup (X t - X T j)l {Ta<oo} 

^ r a <t<oo ' 


= E. r 


— ( X Ta -|- sup X s — ~E x g (X Ta ) 

' 0<s<oo ' 
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We proceed now for r > 0. In this case, we have 

E x G(M)l {M > a} = ^i x G( SUp l{r a <e(r)} 

\ n<c + ^o( ' 7 


0<i<e(r 

= E x G(x Ta + sup ( X t — X ! - a )'j l{ Ta<e ( r )} 

' T a <t<e(r ) ' 

oo 

= E r ! G(X t „+ sup (X Ta+s - X Ta )\re~ rt dtl{ Ta<00 } 

0 <S<t—T a 


— E. T e 


r(X Ta 

[ g(X T a + sup (X Ta+s - X Ta ))re~ rv dvl {Ta<oo} 

Jo ' 0<s<t; ' 


= E,E 0 e- rr ' 
= E x e _rT| 

concluding the proof. 


G(A' Ta + sup Xs)re rv dv l{ Ta<0 o} 
Jo ^ 0<s<v ' 

“E 0 G(X Tq + M) l {Ta<oo} = E x e~ rT ‘ 


a ff(^T a ) l{r a <Oo} J 


□ 


Remark 2. Fluctuation identities as the one presented in the previous Lemma 
in case of exponential or related to exponential functions have been obtained by 


Darling, Liggett and Taylor (1972) for random walks and by Alili and Kyprianou 


(2005) for Levy processes. In case of power functions with positive integer ex¬ 


ponent Novikov and Shiryaev (2004) introduced the Appel polynomials to obtain 


similar identities for random walks, and Kyprianou and Surya (2005) obtained 


the corresponding result for Levy processes. The case of power functions with 
real positive exponent was considered in Novikov and Shiryaev (2007) for both 
random walks and Levy processes. The identity for general functions was ob¬ 
tained by Surya (2007a), see also Surya (2007). 

Lemma 2. Consider a non-negative non-decreasing function f(x ) and a real 
r > 0. Then: (a) The function h{x) = E x f(M) (x G R) is r-excessive, and, 
in consequence, (b) the process {e~ rt h(X t )} is a supermartingale. 


Proof. The fact that (b) follows from (a) is standard, see for example Shiryaev 


(2008). The statement (a) is a corollary of Lemma 2.2 in Christensen et al. 
(2013), as for non-decreasing / we have 


sup f(X t ) = / sup X t = f(M), 

0<t<e(r) \0<£<e(r) J 

concluding the proof. □ 

Proof of the Theore m [1} We finally observe that ([8]) follows from Lemma[l]with 
a = x*. To prove (|9| we observe that V(x) is excessive based on Lemma [2] 
applied to the non-decreasing function G(x) in so, for any stopping time 
r € M, we have 


V(x) > E x e~ rT V(X T ) > E x e~ rT g(X T ). 
This concludes the proof of Theorem [I] 


□ 
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Remark 3. If the equality in ([T]) holds for some x < x*, then, defining the set 

S = {x £ R: V(x) = g(x)j, 


the stopping time 

r** =inf{t >0: X t £5} 

is also an optimal stopping time for the problem § In fact, from the super¬ 
martingale property, as t** < t* , we have 

V(x) > E x e~ rT '*V(X T **) > E x e~ rT *V(X T .), 

obtaining that E x e~ rT V(X T »») = E x e~ rT g(X T *). 

Remark 4. A method to find G*(x ) and x* consists in first imposing condition 
Q for all x G R, i.e. in finding the averaging function of g and M, and then 
finding its largest root. This determines G* for x > x*, in case it is a non¬ 
decreasing function on this half-line. 


The following result gives a sufficient condition in order to verify condition 

Hi- 


Corollary 1 (Surya (2007a)). Assume that there exists a function Q: 
such that 


E X Q(M) = g{x) for all x, 
and a real constant x* such that whenever 


x < x* < y < z 

we have 

Q{x ) < Q{x*) = 0 < Q(y) < Q(z). (11) 

Then G*(x) = Q{x) when x > x* verifies the conditions of Theorem [I] , and 
lim^oo g(x) = 0. 

Proof of the Corollary. Let us check first that 

lim g(x) = 0. (12) 

x—>—oo 


In fact, if limsupt^.Qo g{x) > 0 there exists a decreasing sequence x n —> —oo 
such that g(x n ) > I > 0. But 


g{x„) = E Q Q(x n + M) 

= E 0 Q(x n + M)1 

{x n -\-M>x*} + Eo Q{x n + AL)1 {x^+MKx*} 
< Eq Q(x n + M)\{ Xn+ M> x *} — > 0 as n —> oo, 


by dominated convergence, as Q(x)l{ x j r M>x»} is decreasing in x by hypothesis, 
giving a contradiction, and concluding (12). 

The rest of the proof is immediate as condition (11) implies condition 0- 
In fact, for G defined in ([5]), we have 

V{x) = E X G(M) > E X Q(M) = g{x), 

concluding 0- and the proof of the Corollary. □ 
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3 Polynomial rewards 


Our payoff function is constructed from a polynomial 

p n (x) = x n + a n - 1 x n ~ 1 H-baix, (13) 


where we assume that x = 0 is a root of p n (x). The payoff is the positive part 
of a polynomial, for positive values of the variable x: 


g{x) = ( p n {x + )) + = 


\Pn{ a 

0, 


when x > 0, 
otherwise. 


Observe that, the problem ([2]) with reward function ag(■ + xq) has solution 
aV{- + xo), taking the first coefficient a n = 1 and x = 0 as the smallest root of 
p n in (131 entails no loss of generality for any polynomial with positive leading 


coefficient and at least one root. 


3.1 The averaging polynomial 

We search for a function P n {x) such that 

E x P n {M) = p n (x),x e K. (14) 


It is not difficult to see that this averaging function can be taken to be a poly¬ 
nomial of order n, 

P n (x)=x n +b n _ ix n_1 H- \-biX + b 0 . (15) 


Assume that the first n moments of M are finite and denote them by gk = 
Eo(Al fe ) (fc = 1 ,...,n). Denote po = 1. With this notation, the l.h.s. in 
equation (14), after changing the order in the sums, reads 


n / n 


^ b k Y, Cix e p k _ e = Y[H x\ 


k -o e—o 


£=0 \k—i 


that equating coefficients of equal degree in (14) gives 


Y hCipk-e = ae, l = n, n - 1,..., 0. 


k=i 


This system of equations can be solved recursively backwards, i.e. 


(16) 


— 1 5 

bn— 1 = ^n— 1 Tilli , 

n 

bi — cig ^ ^ bkCfcfik—g) £ — n 2,..., 0, 
k=t+1 


where we put «o = 0. 





Remark 5. An equivalent way to obtain the averaging function P n {x) in (141 
is to write it as 


Pn(x) — ^ ' OkQkix) 


k =1 


where the Qk(x) are the Appell polynomials of the random variable M, intro- 


duced 

in Novikov and Shiryaev (2004 

), applied also in 

Kyprianou and Surya 

(2005 

), Salminen \ \2001 

) and 

Mishura and Tomashyk (2011). 


We have the following simple result. 


Proposition 1. Consider a polynomial p n (x) as in (13). 

(a) The averaging polynomial P n {x) constructed as in (15) has at least one 
positive root. 

(b) If x* denotes the largest root of P n {x), we have p n (x) > 0 for x > x* and 
p n (x) > 0 for x > x*. 


Proof, (a) If P n {x ) has no positive root for x > 0 then P n (x) > 0 for all x > 0. 
As Po (M > 0) > 0, this gives 0 < E 0 P n (M) = p n ( 0) = 0, a contradiction. 

(b) is a consequence of 


p n (x) = E x P n (M) > 0 for x > x*, (17) 

as P X {M > a;*} > 0 and P n (x) > 0 for x > x*. The condition P X (M > 
x) > 0 and inequality P n (x) > 0 for x > x* gives the strict inequality in 
concluding the proof. □ 

Theorem 2. Let p n [x) be a polynomial of degree n with leading coefficient 
a n = 1 and p n { 0) = 0. Define as before 

g(x) = (Pn{x + )) + . 

Denote by P n {x) the averaging polynomial of p n [x) for the random variable M. 

Denote by x* the largest positive root of P n (x). Define G(x) = P„.(x) 

and 

V(x) = F, x G(M), t* = M{t>0: X t >x*}. 

If G(x) is non-decreasing and V(x) > g{x) for x < x*, then, the pair V(x), t* 
is a solution of the optimal stopping problem ([2|. 

Proof. The result follows directly from the application of Theorem |T] □ 


4 Examples 

In order to illustrate our results we first assume that X is the Brownian motion 
and r = 1/2. In this case M has exponential distribution with parameter one (in 
the general case with parameter l/y2r). Its moments satisfy p n = T(n + 1) = 
n\ Observe that for spectrally negative Levy processes, the random variable 
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M is also exponentially distributed, with parameter <I>(r), where <l>(r) is the 
unique positive root of equation ip(\) = r and tjj( A) is the Laplace exponent of 
corresponding Levy process (see for instance Bertoin (1996)). In this case we 
can produce similar results. 


4.1 Example 1: Quadratic polynomials 

Consider p 2 (x) = x 2 + ax. Solving ( jToj ) we obtain 

P 2 (x) = x 2 + (a — 2hi)x + 2(hi) 2 — Hi ~ a/x i 
that has its largest root 


x* = Hi - a/2 + fi 2 - Mi + a 2 /4 = E 0 M - a/2 + y/var 0 M + a 2 /4, 

that is evidently positive that can be checked independently of Proposition [T| 
In case a = 0 we obtain 


x* = Hi — a/2 + hi — Mi + « 2 / 4 = E 0 M + a/ var 0 M 

that gives the solution found in Novikov and Shiryaev (2004). In our particular 
case Hi = 1, Hi = 2, Pi{x) = x 2 + (a — 2)x — a, and x* = 1 — a/2 + ^1 + a 2 / 4. 
For any a £ R. it is evident that G{x) increases after x* and it is not difficult to 
calculate V(x): 


V(x) = {x 2 + ax)l {x>x , } + ((x*) 2 + ax*)e x 

In order to apply Theorem [2| we need only to check the condition V(x) > 
g(x ) for x < x*, but in fact it is only necessary to check this for —a < x < x*. 
Consider the case a < 0, the opposite case is considered similarly. So, we need 
to check the condition 


(x 2 + ax)e x < ((x*) 2 + ax*)e x 


(18) 


for — a < x < x*. The latter inequality holds for x = — a where we have 
the strict inequality and for x = x* where we have the equality. Furthermore, 
function /(x) = (x 2 + ax)e~ x has the derivative f'(x) = —P 2 (x)e~ x which is 
positive between the roots of P 2 (x), the biggest is x*. Moreover, the smallest 
root of P 2 (x) equals 1 — a/2 — sj 1 + a 2 /4 < —a for negative a. It means that 
/(x) increases on (— a,x*) whence we get (18). So, according to Theorem [2j 
(F(x),x*) create a solution of ([2|. The same is true for a > 0. In Fig. |4.1| we 
plot the solution for a = — 1 and a = 1. 


4.2 Example 2: Cubic polynomials 

Consider ps(x) = x 3 + ax 2 + bx. Solving © we obtain 

P 3 (x) = x 3 + (a — 2>hi)x 2 + (b — 3 /x 2 — 2(a — 3 hi)hi) x 

- H 3 - (a - 3 hi)hi ~ (b - 2 >H 2 - 2(a - 3hi)hi)^i- 
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Figure 1: Example 1 with a = 1 (left) and a = — 1 (right). Here p 2 (x) is plotted 
continuously, P 2 (x) (dashed) gives the roots. The thick lines are the respective 
solutions V{x). Observe that in case a = 1 (left) the averaging function P 2 does 
not remain non-positive for values smaller than the root x* ~ 2.62. 


If we further assume that X is the Brownian motion and r = 1/2, we have 
P 3 (x) = x 3 + (a — 3)x 2 + (b — 2 a)x — b 

and 

V(x) = (x 3 + ax 2 + bx)l{ x>x .y + ((x*) 3 + a(x*) 2 + bx*)e x ~ x l{ a < x .}. 

If b > 0, P 3 evidently has at least one positive root since P 3 (0) < 0 and 
P 3 (+oo) = + 00 . Let b < 0, a > —2, then P 3 ( 1) = — a — 2 < 0, and at least one 
positive root exceeds 1. If b < 0, a < —2, then P(—a) = —a 2 — 6(1 + a) < 0 
and at least one positive root exceeds —a. So, in any case P 3 has positive 
roots, in accordance with Proposition [l] but we have checked this indepen¬ 
dently. Now, in order to apply Theorem [2j consider some particular cases. 
In the case when 3 < a < | A 8, (the case of positive coefficients, for exam¬ 
ple, a = 4, b = 10) p 3 (x) has only one root x = 0 because other roots that 

should equal Xi i2 = — | ± \J^ — b do not exist (discriminant is negative, since 

N— ^ ^ - 6 < 0), the derivative P/(x) = 3x 2 + 2 (a — 3)x + b — 2a has two 

negative roots therefore it is positive on [0, 00 ), and P 3 (x) increases on [x*,oo), 
even more, it increases on [0, 00 ) being negative on [0,x*). Moreover, as in the 
example with quadratic polynomials, we need to check inequality 

(x 3 + ax 2 + bx)e~ x < ((x*) 3 + a(x*) 2 + bx*)e~ x 

on the interval [0, x*] but on the interval [0,x*] the derivative of the function 
(x 3 + ax 2 + bx)e~ x being equal —P 3 (x)e~ x is positive therefore both conditions 
of Theorem [ 2 ] hold. In the case a = b = 0, i.e. p 3 (x) = x 3 , we have 

P 3 (x) = x 3 — 3x 2 , 

with largest root x* = 3 (see Fig. [ 3 ]). Evidently, P 3 increases on [x*, 00 ) because 
its derivative 3x 2 — 6x is positive on the interval [2, 00 ) and —P 3 is positive on 
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Figure 2: Example 2 with a = b = 0 (left) and a = —9/8, b = 3/8 (left). Here 
P 3 {x) is plotted continuously, Q 3 (x) (dashed) gives the root x* in each case. 
The thick lines are the respective solutions V(x). 


(0,3) which supplies both conditions of Theorem [2j An example for Brownian 
motion with r = 1/2 and polynomial with positive b and negative a is shown in 
Fig- @ We put in this case 

p 3 (x) = x 3 - (9/8)x 2 + 3/8, P 3 (x) = x 3 - y x 2 + y x - 

and the largest root is x* = 3.3815. 


4.3 Example 3: Kou’s process for a quartic polynomial 


A diffusion process {X t } with two sided exponential jumps, defined by the 
formula 

N t N' 

x t =at + aW t + J2 Y k-J2 Y k’ 


k=1 


fc =1 


is known in the financial literature as a Kou’s process (see Kou (2002) and Cont 


and Tankov (2004)). Here {Nt} (resp. is a Poisson process with param¬ 


eter p (resp. v) and {Yfc} (resp {Y/}) is a sequence of independent exponential 
random variables with parameter a (resp /?). The cliaractristic exponent ([Tj) of 
the process is given by 

i/j{z) =az+ \cf 2 z 2 + p— - 

2 a — z z + p 

and the density of the maximum M in this case is a mixture of two exponentials 
Jm{x) = Airie~ riX + A 2 r 2 e~ r2X , 


with coefficients 


1 ~ n//3 = 1 - r 2 /P 

1 — ri/r 2 ’ 2 l-r 2 /ri’ 
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where 0 < n < r 2 are the positive roots of the equation if(z) = r (see Mordecki 
(2003)). In consequence the moments are given by 


Mfc = k\ -f-+ 


A i A, 


We consider a quartic polynomial 


If we denote P(x) 


p(x) = x 4 + a^x 3 + a 2 x 2 + a±x. 
x 4 + 63 a ; 3 + b 2 X 2 + b\x + bo, applying (16) we obtain 


h = a 3 - 4/ii 

62 = a 2 - 6/x 2 - 363/Lti 

61 = ai — 4 /i 3 — 363^2 — 2&2/T1 
= — (/X4 + &3/X3 + 6 2 /i2 + &1/T1) 


Assuming that there exists a value a;* that satisfies the conditions of Theorem 
[2j we write the possible value function 

V(x) = E x P(x + M) 1 {X+M > X . } = B ie r ^ x ~^ + B 2 e r ^ x - x '\ 

where 

pOO pOO 

Bi = Airi / P(z + x*)e~ riZ dz, B 2 = A 2 r 2 / P(z + x*)e~ r2Z dz. 

Jo Jo 

To proceed we choose values for the parameters: 


a = 2 , cr=l, n = 1 , 1 / = 1 , a = 2 , /? = 2 , r = 6 . 
and choose the polynomial 

p(x) = a: (a; — l)(x — 2)(x — 3) = x 4 — 6 x 3 + llx 2 — 6 x. 

We obtain n = 1.4327, r 2 = 2.8740, giving A± = 0.5656, A 2 = 0.4344. In 
consequence x* = 4.3706. For this sets of parameters the conditions of Theorem 
[ 2 ] are fulfilled (see Figure [3]). 
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